Planar Skyrmions at High and Low Density 



R. S. Ward* 

Department of Mathematical Sciences, 
University of Durham, 
Durham DH1 3LE 

February 1, 2008 



Abstract 

The 0(3) Skyrme system in two space dimensions admits topological soliton solu- 
tions. This paper studies the transition between the high-density crystalline phase of 
such solitons and the low-density phase where there are multi-Skyrmions localized in 
space. The details depend crucially on the choice of the potential function. Two such 
choices are investigated: in the first system, multi-Skyrmions at low density form a 
ring; while in the second, there is an explicit crystal solution, and the preferred low- 
density configurations are chunks of this crystal. The second system is a particularly 
good analogue of three-dimensional Skyrmions. 
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1 Introduction 



Topological solitons are of interest both mathematically and in many areas of physics 
(as models for particles, topological defects in condensed-matter physics, etc). A par- 
ticularly important question is the nature of, and transition between, the high-density 
and low-density phases of such solitons. Ultimately, one would like to understand the 
thermodynamics of these systems; the starting-point for this is to investigate their 
zero-temperature behaviour. At its most basic, this amounts to considering a fixed 
number N of solitons confined to a finite volume, and investigating the way in which 
static classical A-soliton solutions depend on the volume (or density). Alternatively, 
one may use the dimensionless ratio (size of soliton) / (size of space) as a parameter. 

The typical situation is as follows. At large density, there is a high degree of 
symmetry and uniformity, and in particular it may not possible to identify individual 
solitons; in Skyrme models, for example, one gets a periodic crystal-like structure of 
half-Skyrmions. At low densities, by contrast, solitons (or multi-solitons) become 
localized in space, and there is less symmetry. The localized multi-solitons may 
be 'chunks' of the high-density crystal, or may have a quite different (for example 
shell-like) shape — this depends on the details of the system. 

The picture for the three-dimensional Skyrme model may be summarized as fol- 
lows. Let N denote the Skyrme number, and E the normalized energy-per-Skyrmion 
(so the Bogomolny-Faddeev bound is E > 1). At high density, the ground state is a 
triply-periodic lattice of half-Skyrmions E3 lU 03 E] , with energy E = 1.038 at its 
most favourable density. At low density (for example for an isolated multi-Skyrmion 
in R 3 ), and for relatively low values of N (in particular for N < 22), the minimal- 
energy static Skyrmions take the form of polyhedral shells [7]. Their normalized 
energy E is a decreasing function of N, and it appears [B] that these polyhedra have 
energy E ~ 1.06 for large N. Since this is larger than the energy of the Skyrme 
crystal, one expects that there is a critical value N c such that for N > N c , the 
minimal-energy Skyrmion resembles a chunk of the Skyrme crystal. Such a lattice 
chunk will have energy E ~ 1.038 + kN^ 1 / 3 , with the second term being a surface 
contribution. There have been attempts [HI to estimate the constant k, and hence to 
determine N c , but these have not been definitive. There are also other possibilities 
for the shape of large- A Skyrmions, for example a multi-shell structure |10j . 

The present paper deals with the two-dimensional 0(3) Skyrme system, which 
one may view as an analogue of the three-dimensional well as being of 

interest in its own right. The 2-D Skyrme system contains, as a limiting case, the 
two-dimensional 0(3) sigma model; the sigma- model crystal, obtained by imposing 
periodic boundary conditions, has been studied both as a model quantum field theory 
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jllj and in connection with the dynamics of (classical) solitons |121 I13| . But sigma- 
model solitons do not have a fixed size, and in particular tend to shrink and decay; 
so in that sense they are not true solitons. We shall consider only the case where 
both the Skyrme term and a potential term V are present, and the soliton size is 
consequently fixed. The interest here is in the nature of the low-density and the 
high-density configurations, and the transition between them; as we shall see, the 
details of these depend crucially on the choice of potential V. Two different systems, 
corresponding to two different choices of V, will be investigated in detail. 



2 Two- Dimensional Skyrmions 

The two-dimensional Skyrme system involves a unit 3- vector field (ft = (0i, 4>2, (fts) 
defined on (2+l)-dimensional space-time. We are interested only in static configu- 
rations, so (ft depends on the spatial variables x 1 = (x l ,x 2 ), thought of as (local) 
coordinates on a two-dimensional space S. If S is a compact surface, then (ft has a 
winding number N, which we think of as the number of Skyrmions. One special case 
is where space is the plane R 2 with boundary condition eft — > (0, 0, 1) as r — > oo - 
this corresponds (by conformal invariance) to the 2-sphere S = S 2 . Another case 
is that of periodic boundary conditions <p(x , x 2 ) = ^(x 1 + L,x 2 ) = ^(x 1 ^ 2 + L), 
which corresponds to the 2-torus S = T 2 . In either of these cases, we may regard 
the spatial metric as being the flat (Euclidean) metric, and we shall do so in what 
follows. 

The integer N can be either positive or negative; without loss of generality, we 
shall assume N to be positive. The energy density of a configuration (ft is 

£ = h( d j& ■ ( d 3$) + l«^ 2 + \olV{$), (1) 

where = (ft ■ d\(p x 82(f), V is some potential function, and a is a dimensionless 
constant. The length-scale in this system is determined by the ratio of the coefficients 
of the f2 2 and V terms, and so is fixed by choosing these coefficients to be equal as 
in (jT]). Static multi-Skyrmion solutions are critical points of the (normalized) energy 
functional 

E=-±- I Sd 2 x. (2) 
4vriV J v ; 

There is a topological (Bogomolny) lower bound on the energy, namely 

E > 1 + j- [ s/v^dw, (3) 

where did is the usual area element on the space S 2 of unit vectors (in other words, 
J du = 4tt). Under certain circumstances, this bound can be saturated; the following 
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statement is adapted from [Tl], see also ^2]- Write z = x 1 + ix 2 , and let W = 
(01 + /(l — ^3) denote the stereographic projection of 0. Let W(z) be a complex- 
analytic function satisfying a first-order ordinary differential equation of the form 
dW/dz = F(W) for some function F. Then the configuration M^(z) saturates the 
bound ©, and hence is a static Skyrmion solution, provided that the potential V is 
given by 

v= ^\f{wT 

(i + |iy| 2 ) 4 K) 

The simplest example of this is where F is a constant; the corresponding system 
(on R 2 ) has been investigated in some detail, for example in |16l I17j , In this case, 
there is a repulsive force between Skyrmions, and consequently there are no static 
multi-Skyrmion solutions. 

For any system of the form Q , the nature and shape of multi-Skyrmion solutions 
depend on the choice of V, and many different choices have been considered |14l 
118] . The following two sections will deal with two possible choices. They are both 
analogous to the three-dimensional Skyrme system, in that they allow both crystal- 
like and ring-like solutions (the latter being the counterpart of the polyhedral shells 
mentioned earlier). 

To conclude this section, it is worth mentioning yet another case which has been 
extensively investigated, namely V(<p) = 1 — 03 ^HEOl- In this case, the static 
multi-Skyrmion solutions on R 2 appear to form a lattice-like structure |19[ I18| — for 
example, for even values of N one gets a lattice of double-Skyrmions as the lowest- 
energy state. There are also other local minima of the energy, but there do not 
appear to be any ring-like solutions (except for N = 2). The thermodynamics of this 
system has been studied numerically |21j . 



3 The System V = 1 - <#j 

In this section, the potential function V appearing in Q is taken to be V((f>) = 1 — 03- 
Let us discuss, first, the localization-delocalization transition for this system. A useful 
order parameter in this regard is the quantity 




where A is the area of the 2-space S. At high density, there is a high degree of 
symmetry, and we expect to have (03) = 0; while at low density, the field will 
localize in space and we will have (03) ^ 0. The transition between these two phases 
occurs at some critical density p c . This transition was investigated by taking space 
to be a sphere (S = S 2 ) in m particular, this paper looked at (03), as a function 
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of A, for a single skyrmion on S 2 . Using approximate analytic expressions for the 
Skyrmion indicated that the transition occurs at A = 47r\/5, ie the critical density is 

p c = l/(47iV5) = 0.036. (5) 

Numerical simulations gave results that supported this estimate. 

Let us now compare this to the situation on the torus T 2 . So we are looking for 
doubly-periodic solutions — periodic in both x 1 and x 2 with period L. If a = (ie 
for the 0(3) sigma-model) , then solutions correspond to elliptic functions, and these 
have topological charge N = 2 in the unit cell. So we expect the minimal-energy 
crystal in the Skyrme case also to have an N = 2 cell. Such solutions can be found by 
numerical minimization of the energy functional. It turns out there are two different 
crystals, namely one with (approximate) symmetry A: 

Al. (x,y) i-» {-x,y), (fa, fa, fa) >-> (fa, -fa, fa), 

A2. (x,y) i ^ (y,x), (fa, fa, fa) ^ (~ fa, fa, fa), 

A3, (x, y) i-> (x + L/2, y), (fa, fa, fa) i-> (0 3 , -0 2 , 0i); 

and one with symmetry £?: 

Bl. (x,y) 1 ^ (-x,y), (fa, fa, fa) ^ (fa,~<l>2,fa), 
B2. (x,y) 1 ^ (y,x), (fa, fa, fa) >-> (~fa,fa,~fa), 
B3. (x,y) 1 ^ (x + L/2,y), (fa, fa, fa) i-> (-fa, -fa, fa)- 

In the a = system, symmetries j4 and -B are equivalent; but this degeneracy is 
broken by the potential term V. The situation is illustrated in Figure ^ which plots 
the normalized energy E for the two solutions, and the quantity \£' = (fa) a for the 
type-^4 solution, as functions of the periodicity L. The value of the parameter a used 
in obtaining these graphs is a = 0.1. The Bogomolny bound (j3J) for this system is 
E > 1 + 7ra/4 = 1.0785. The lowest crystal energy is attained when L = Lq = 6.0, ie 
at density p = 0.056; this solution has symmetry of type B, and normalized energy 
-Ecrys = 1-088. The type-i? solution exists for lower density, and is always delocalized, 
ie has (fa) = 0. 

At relatively high density, in fact for L < L\ = 7.0, the solution with type-^4 
symmetry also has (fa) = 0; but it is unstable, and its energy (which is not plotted 
in Figure ^) is higher than that of the type-i? solution. When perturbed, it decays 
to the type--B solution. For L > L\, the type- A solution delocalizes (^ 7^ 0) and 
becomes stable. (It no longer has the symmetry A — in fact A3 is invalid — but we 
still refer to it as the type- A solution.) This is illustrated in the right-hand diagram 
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Figure 1: E and \l/ = (fa) a for the periodic N = 2 solutions. 



of Figure Q the numerical results indicate that 

lim * Kt 0.06, 

but this limit is difficult to study owing to the instability which sets in below L\. In 
the range L\ < L < L<i = 7.35, the energy of the type- A solution is greater than that 
of type B — in other words, there are two local mimina of the energy, one of which 
(type B) is more symmetric than the other. For L > L2, the minimal-energy N = 2 
solution is the localized type- A one. 

We see, therefore, that there is a transition between a dense, highly-symmetric 
homogeneous phase and a low-density, less homogeneous phase. The latter appears 
at L\ = 7.0, corresponding to a density p = 2/7 2 = 0.041, and it becomes the 
minimal-energy state at L2 = 7.35, corresponding to the critical density 

p c = 2/7.35 2 = 0.037. (6) 

The close agreement between © (for N = 2 Skyrmions on T 2 ) and (J5J) (for N = 1 
Skyrmions on S 2 ) is remarkable. 

Let us now consider the zero-density limit L — > 00. The only known static 
solutions on R 2 are rotationally-symmetric rings [22] . The energy of these goes like 
E ~ j3 + 7/iV 2 , where (3 and 7 are constants. The significant fact here is that j3 
is lower than E crys ; this has been checked numerically for a < 1. For example, for 
a = 0.1, we get (3 = 1.082, compared with -E crys = 1.088. The implication of this is 
that, unlike the case for three-dimensional Skyrmions, isolated crystal chunks are not 
energetically favourable, even for large values of N. In fact, it seems likely that the 
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minimum-energy configuration on R 2 is, for all values of N, a single 0(2)-symmetric 
ring (other local minima may, however, exist). In this respect, the system is different 
from the three-dimensional Skyrme model. 

4 The System V = (1 - 0§)(1 - <fi) 

The motivation for the potential V studied in this section comes from requiring 
that the Bogomolny bound © be saturated by a doubly-periodic solution. As a 
consequence of this, the energy of the corresponding Skyrme crystal will be as low as 
it possibly can be; and so an isolated multi-Skyrmion will, for N large enough, take 
the form of a crystal chunk rather than a ring-like configuration. In this respect, it 
will be a better analogue of the three-dimensional Skyrme system. 

Our assumption, therefore, is that the system should admit a solution W(z) which 
is an elliptic function. One particularly simple choice is to take W(z) = 2p(z), where 
p is the Weierstrass p- function with parameters g2 = 1 and 53 = 0. So the function F 
appearing in @ is given by F(W) 2 = 2W(W 2 — 1); and the corresponding potential 
function is easily seen to be 

F(^) = 16(l-0 2 )(l-</> 2 ). (7) 

For the rest of the section we shall use this V . 

The fundamental cell of the lattice is a square with sides of length Lq = 2 K (0.5) = 
3.708, where K(m) is the complete elliptic integral of the first kind with parameter 
m. This cell contains two units of topological charge (that is, the map from T 2 
to S 2 has winding number N = 2); the distribution of charge- and energy-density 
indicates that it should be thought of as a square lattice of half-Skyrmions. It has 
the symmetry of type A defined in the previous section. If we scale the solution to 
fit a lattice with edge- length L, then we get a solution with normalized energy 

E(L) = 1 + 2.3729 X | (Jj- + . (8) 

If L = Lq, the solution saturates the lower bound lj3*j). which in this case is E > 
1 + 2.3729 a. 

For large L (low density) one expects that this solution (a scaled Weierstrass 
function) will no longer be the minimal-energy doubly-periodic N = 2 solution. In 
fact, the minimal solution will correspond to an isolated N = 2 lump at some point 
on T 2 — and so will have (^3) 7^ 0. The result of a numerical investigation of this 
phenomenon is presented in Figure I^J it was obtained by finding local minima of 
the energy with N = 2 for a range of values of the periodicity L. The value of 



7 



1.35 



1.34 
1.33 
1.32 



w 1.31 
1.3 
1.29 
1.28 
1.27 

5 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 6 

L 

Figure 2: Energy E for the two N = 2 solutions. 

a was taken to be a = 0.1. For L < 5.15, there is only one solution, namely the 
Weierstrass-type one with (^3) = 0. Its energy grows fairly rapidly with L, but it 
remains a stable solution (ie a local minimum), at least up to L = 20. For L > 5.15, 
another solution (local minimum) appears — a spatially-localized one with (^3) 7^ 0. 
Its energy grows less rapidly with L, and for L > 5.73 it is the lowest-energy solution. 
In this sense, there is a phase transition at L\ = 5.73, corresponding to the critical 
density p c = 2/5. 73 2 = 0.061. It would be interesting to study this transition in 
other ways, for example by examining the system on a sphere (as in |231 115j ) or by 
Monte-Carlo simulations at finite temperature (as in |21j). 

Let us now look at the L —* 00 limit, and consider finite- iV solutions on the 
infinite plane R 2 , with the boundary condition ^3 = 1 at spatial infinity. First, the 
system admits ring-like solutions. For these, we have ^3 = —1 at a single point O, 
and 03 = on a deformed ring centred at O. Because of the factor (1 — 4>\ ) in 0, this 
ring is not quite rotationally-symmetric about O; its symmetry group is a subgroup 
of the rotation group 0(2), namely the dihedral group D2N- In other words, these 
solutions resemble polygonal rings, analogous to the polygonal shells of the three- 
dimensional Skyrme system. A numerical investigation (for the case a = 0.1) reveals 
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that their normalized energy E is very well fitted by the function 

E = 1.262 + 0.227/7V 2 (9) 

for N > 2. For example, the energy of the N = 8 ring is E = 1.266; for comparison, 
the topological lower bound is E = 1.2373. This N = 8 ring is depicted in the left- 
hand plot of Figure it is a stable static solution with Die symmetry, and is best 
thought of as a ring of sixteen half-Skyrmions. 




Figure 3: Contour plots of 02 for two static N = 8 solutions of the system (j7J). 

The other localized solution illustrated in Figure El corresponds to a clump of the 
Skyrme crystal — in the right-hand plot, we see a 2 x 2 clump of crystal, again con- 
taining sixteen half-Skyrmions. The energy of such a clump is approximately equal 
to the (bulk) energy of the crystal plus a surface contribution. If the circumference 
of the clump has length MLq, where M is a positive integer, then one expects that 
the normalized energy will have the form 

EttE + M5E/N (10) 

where Eq = 1 + 2.373 a and where 5E is the contribution from a single edge of lattice- 
cell. One may estimate 8E as follows (this is analogous to the discussion in J]). Start 
with a single cell, with coordinate range 0<x<Lo,0<y<Lo- Now stretch the cell 
in the positive x-direction, so that the x-range becomes < x < do. The boundary 
conditions are that the field (ft remains periodic in the y-direction, maintains its 
original lattice-configuration on the edge x = 0, and tends to its asymptotic value 
(0,0,1) as x — » oo. We can then relax this configuration, minimizing its energy 
numerically to obtain a mimimum energy E\. The deformed cell still contains two 
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units of Skyrme charge, so the resulting estimate for 6E is 

SE ps E x - 2E . (11) 

For example, if a = 0.1, one gets 5E = 0.026. Consequently, our estimate for the 
energy of a localized square lattice clump of charge N, with a = 0.1, is 

E ps E + 2^2 5E/VN = 1.237 + 0.074/ VN; (12) 

but for low N this will be an underestimate, since there are also corner effects to be 
considered. These corner effects presumably give an additional contribution to E of 
the form C/N, where C is a constant. 

If we take a single lattice cell (N = 2), then it resembles a square array of four 
half-Skyrmions, but this amounts to the same thing as a ring of four half-Skyrmions 
— indeed, the N = 2 solution on R 2 has a unique shape. For N = 8, however, there 
is the possibility of the two solutions depicted in Figure 01 and each of these is in fact 
realized, as a local minimum of the energy. The formulas (|12ft and © suggest that 
the energy of the N = 8 ring should be higher than that of the square clump, by an 
amount AE = 0.0027. The numerical result is that AE = —0.0035, which gives an 
estimate for the corner contribution which has to be added to (|12|) . namely 0.05/iV. 

For larger crystal chunks, the ring will be less favourable. For example, we expect, 
from Q and the modified version of Q12j) . that the minimal-energy N = 18 solution 
will correspond to a 6 x 6 square array of half-Skyrmions; and in particular, that the 
energy of the N = 18 ring will be higher than that of the 6x6 clump by an amount 
AE = 0.0057. It should be emphasized, however, that there are likely to be many 
local minima, of which the square chunk and the ring are just two examples. All this 
is very similar to the situation for the three-dimensional Skyrme system. 

In conclusion, this particular planar system, in addition to admitting an explicit 
crystal solution which saturates the topological lower bound on the energy, is a good 
model for three-dimensional Skyrmions. 
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